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ABSTRACT 

Kraichnan (1965) proposed that MHD turbulence occurs as a result of 
collisions between oppositely directed Alfven wave packets. Recent work has 
generated some controversy over the nature of non linear couplings between 
colliding Alfven waves. We find that the resolution to much of the confusion lies 
in the existence of a new type of turbulence, intermediate turbulence, in which 
the cascade of energy in the inertial range exhibits properties intermediate 
between those of weak and strong turbulent cascades. Some properties of 
intermediate MHD turbulence are: (i) in common with weak turbulent cascades, 
wave packets belonging to the inertial range are long lived; (ii) however, 
components of the strain tensor are so large that, similar to the situation in 
strong turbulence, perturbation theory is not applicable; (iii) the breakdown of 
perturbation theory results from the divergence of neighboring field lines due to 
wave packets whose perturbations in velocity and magnetic fields are localized, 
but whose perturbations in displacement are not; (iv) 3-wave interactions 
dominate individual collisions between wave packets, but interactions of all 
orders n > 3 make comparable contributions to the intermediate turbulent 
energy cascade; (v) successive collisions are correlated since wave packets are 
distorted as they follow diverging field lines; (vi) in common with the weak 
MHD cascade, there is no parallel cascade of energy, and the cascade to small 
perpendicular scales strengthens as it reaches higher wave numbers; (vii) For an 
appropriate weak excitation, there is a natural progression from a weak, through 
an intermediate, to a strong cascade. 



Subject headings: ISM: general - MHD -turbulence 
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1. 



INTRODUCTION 



There appears to be some consensus regarding the notion that turbulence in the ionized 
interstellar medium is of magnetohydrodynamic (MHD) origin. However, even three decades 
after Iroshnikov (1963) and Kraichnan (1965) first, independently, presented their ideas on 
MHD turbulence (in an incompressible, highly conducting fluid), there is much debate on 
what this theory really is. Recent controversy centers on the existence, or otherwise, of 
certain nonlinear interactions among Alfven waves. Sridhar & Goldreich (1994; hereafter 
SG) have argued that there are no resonant 3-wave interactions in weak MHD turbulence. 
They also asserted that the Iroshnikov-Kraichnan (IK) theory is incorrect, and constructed 
a theory of weak MHD turbulence based on resonant 4-wave interactions; in this theory, 
non linear interactions strengthen on small spatial scales, resulting ultimately in a strong 
cascade proposed by Goldreich & Sridhar (1995; hereafter GS). However, Montgomery 
& Matthaeus (1995) have maintained that resonant 3-wave interactions are non empty, 
holding these responsible for the anisotropic cascade seen in the numerical simulations 
of Shebalin, Matthaeus, & Montgomery (1983). Moreover, Ng & Bhattacharjee (1996; 
hereafter NB) have recently presented calculations which show that small amplitude wave 
packets do interact via 3-waves, so long as the = {z is the direction of the local, mean 
magnetic field) Fourier components of velocity and magnetic field perturbations are non 
zero. The present investigation is an analysis and resolution of this controversy. 

We consider magnetic turbulence in an incompressible fluid, although we simply call it 
MHD turbulence. When the density and transport properties of the fluid are constant, the 
equations of (incompressible) MHD read 



where v is the velocity, b = 'B/^/A^Tp is the magnetic fleld in velocity units, and p is the 
ratio of total pressure to the density. The dissipation provided by k and 7 is important 
only on small spatial scales, so we ignore the dissipative terms. The homogeneous state, 
vq = 0, Bo = Bqz, is a stable, static solution of equations (@). Shear Alfven waves and 
pseudo Alfven waves are the two kinds of linear perturbations about this equilibrium, the 
latter being the incompressible limit of the slow magnetosonic wave. Both kinds of waves 
have the same dispersion relation, namely, u = v^lkzl, where va = Bq/ ^/A^Tf) is called the 
Alfven speed. The perturbed velocity and magnetic flelds are related by 5w = ±6h, where 
the upper/lower signs correspond to waves traveling antiparallel/parallel to Bq (with kz < 



dth 



Vx(v X b) + /tV^b, 

-(v ■ V)v + (b ■ V)b - + 7VV, 

V-v = V-b = 0, 



(1) 
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and kz > respectively). Equations (0), with k = 7 = 0, possess the remarkable property 
of allowing for nonlinear generalizations of the linear Alfven waves. Mutual cancellation of 
nonlinear terms permits the following wide class of exact solutions: if 5v(x) = — 5b(x) at 
some instant of time, t = 0, it can be checked that 6\'{x,y, z — VAt) = —6h{x,y, z — VAt) 
for all time, irrespective of the functional form of (5v(x) (c.f. Parker 1979). This nonlinear 
solution describes a wave packet of arbitrary form traveling nondispersively in the direction 
of Bq. Similarly, we can also construct another class of nonlinear solutions, with 6v = 6h, 
that travels nondispersively in a direction opposite to Bq. Both types of nonlinear solutions 
are stable, and the dynamics is simple so long as there is no spatial overlap ("collisions") 
between oppositely moving wave packets. 

Let us consider a situation that is less restrictive than perturbations about a static 
equilibrium (a uniform magnetic field). We imagine that turbulent motions on a large 
scale are set up in the fluid by stirring it in a random but statistically steady fashion. 
These create large-scale, disordered velocity and magnetic fields. Kinetic and magnetic 
energies will cascade to smaller spatial scales.Q The central problem of MHD turbulence 
is to determine the statistical steady-state amplitudes of the fluctuations in v and b, on 
intermediate spatial scales — the so called inertial-range power spectrum. Mathematically 
speaking, this cascade through spatial scales should emerge from the effect of the nonlinear 
terms in equations (^. Physically, it is common to speak of interactions between "eddies" 
(c.f. Frisch 1995). Do interactions between eddies of dissimilar sizes make significant 
contributions to the form of the cascade? For hydrodynamic turbulence, the answer seems 
to be, "no, the dominant interactions are between eddies of similar spatial scales". The 
reason for such a locality in interactions between eddies is that the sweeping due to the 
velocity field of a large sized eddy (on a smaller eddy) may be transformed away by a local 
Galilean transformation. Kraichnan noted that, for MHD turbulence, such a transformation 
has no effect at all on the magnetic Geld of the large eddy; the magnetic field of a large 
eddy acts upon smaller eddies in much the same manner as the mean magnetic Geld does 
on Alfven wave packets. Hence MHD turbulence in an incompressible Guid should reduce 
to turbulence of interacting Alfven wave packets. We recall from the previous paragraph 
that finite amplitude wave packets that travel in the positive/negative 2;-directions do so 
without change in form, so long as oppositely directed packets do not overlap. Kraichnan 
realized that the cascade of energy in MHD turbulence occurs as a result of collisions 
between oppositely directed Alfven wave packets. 

To derive the results of the IK theory, consider a statistically steady, isotropic excitation 



^On scales small enough for the wV^b and the 7V^v terms to be important, kinetic and magnetic energies 
will dissipate into heat. 
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of amplitude <^ va, on outer scale £, of the static equilibrium mentioned earlier. Alfven 

waves have Sv ~ Sb, so that vg is the amplitude of excitation of both velocity and magnetic 
fields. The resulting fluctuations may, at any time, be decomposed into Alfven wave packets 
with scales \ ^ i traveling in the positive and negative z directions. IK assumed that the 
energy transfer is local and isotropic in k-space. Collisions between oppositely directed 
wave packets occur over times of order cu^^ ~ {vAk)~^, and these create small distortions. 
During one collision, each wave packet suffers a fractional perturbation 



During successive collisions, these perturbations add with random phases. The number of 
collisions for the fractional perturbations to build up to order unity is 



The energy cascade rate is e ~ vf/tx, where the cascade time is given by 

~ Nx/{vAk) ~ VA/{kv\). The three-dimensional energy spectrum, E{k), is 
related to the velocity fluctuation by v\ ~ k^E{k). Using Kolmogorov's hypothesis of the 
scale independence of e, we obtain the following scalings for the inertial-range of the IK 
theory: 




(2) 




(3) 




m 




(4) 



With these scalings. 




(5) 



so the cascade weakens as A decreases. 



2. 



INTERMEDIATE TURBULENCE 



2.1. 



A New Cascade Based on 3- Wave Interactions 



SG argued that the IK theory, although seemingly plausible, was basically incorrect. 
Here we discuss the reason for the failure of this theory, and propose a new cascade based 
on 3-wave interactions. Let us begin by listing three key features of the IK theory. 
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(i) Wave packets of size A live for A^^ wave periods; large values of Nx correspond to weak 
interactions between wave packets of size A. From equations (0) and (0), it may be verified 
that A^A oc 1/A^/^, so that the cascade weakens as it progresses into the inert ial-range. 

(ii) From equation (H), the fractional perturbation suffered by a wave packet during one 
collision is ~ {vx/va) ^ 1; In the IK theory, interactions between wave packets are 
described by the lowest order non linear terms. 

(iii) Isotropy is assumed in the derivation of the IK theory. 

The derivation of this theory is essentially heuristic, roughly as given above in equations (H) — 
(1). We note this, so that it is clear to the reader that there isn't a more rigorous version of 
the "theory" that we happened not to mention. Together, features (i) and (ii) imply that 
the "elementary interactions" between Alfven waves must satisfy the 3-wave resonance 
conditions: 



where = VA\kz \ ■ Shebalin, Montgomery & Matthaeus (1983) noted that the only non 
trivial solutions of equations (||) require that one of either kiz, or must be zero. This 
implies that waves with values of kz not present initially cannot be created during collisions 
between oppositely directed wave packets. As they point out, there is no parallel (i.e. along 
kz) cascade of energy. Thus the turbulence must be anisotropic, and energy should cascade 
to large k±. 

In what follows, we derive an anisotropic version of IK theory. As before, we imagine 
that the system is stirred in a statistically steady and isotropic fashion such that vi -C va 
on outer scale i. The absence of a parallel cascade imphes that wave packets belonging to 
the inertial range have parallel scales i and perpendicular scales A -C £. We estimate the 
spectrum of the anisotropic cascade by modifying the arguments of equations (0) — (|) so 
as to keep track of parallel, and perpendicular scales. In one collision between two such 
oppositely directed wave packets, the fractional perturbation is given by 



ki + k2 = ks , 



(6) 



5vx 



dvx 
dt 



{vAkzVx) 



1 



^Vx 



< 1. 



(7) 



\va 



Adding up the perturbations due to successive collisions with random phases, the 
number of collisions over which the perturbations grow to order unity is 
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Making use of the steady-state relation ( "Kolmogorov's hypothesis" ) 

-e, (9) 



we obtain scahng relations for the velocity fluctuations, as well as the three dimensional 
energy spectrum of the anisotropic, 3-wave cascade 



E{KM)^^- (10) 



It follows that 



--(^n^ (n, 

so nonlinearity increases along the inertial range of the new cascade. 

The anisotropic cascade differs from the original IK cascade by: (i) being anisotropic; 
(ii) having a different spectrum; (iii) strengthening at high wave number. This final 
difference has a profound consequence. It implies that the anisotropic cascade has a limited 
inertial range, thereby diminishing its applicability in astronomical contexts where the 
excitation at the outer scale is likely to be quite strong. 

It turns out that this anisotropic cascade is an example of a new type of turbulence, 
which we call Intermediate Turbulence, because it has properties intermediate between 
those of weak and strong turbulence. In particular, intermediate turbulence does not submit 
to perturbation theory. 



2.2. The Failure of Perturbation Theory in Intermediate Turbulence 

2.2.1. Field Line Geometry 

To lowest order in perturbation theory, wave packets move along field lines. Thus the 
breakdown of perturbation theory may be understood physically by studying the geometry 
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of the divergence of a bundle of field lines. Assume that the mean field lies along the z 
axis. Consider wave packets localized in velocity and magnetic field perturbations, but not 
in displacement, having longitudinal scale i and transverse scale A with Xji < 1. Let us 
require that the velocity fluctuation, vx, is small enough, 

X- 1^ « 1. (12) 

so that Nx <^ 1 (cf. eqn. P]). The rms differential inclination of the local field across scale 
A is 6'a ~ vx/va- It is correlated over distances I and A parallel and orthogonal to z. 

Let us focus on a pair of neighboring field lines separated by A at z = 0. The 
expectation value of the separation between these field lines. A, varies such that 



A2~A2 + 02^1^1, (13) 

for \z\ ^ The distance along z, over which |A| increases by a factor of order unity from 
its initial value of A at 2; = 0, is 

is an function of A. The significance of L^, in intermediate turbulence follows because 
turbulence involves the transfer of energy across scales. If x 1, single interactions 
between wave packets result in small perturbations. Cascading of energy requires of order 

iVA ~ ^ ~ ^ » 1 (15) 

such interactions. 



2.2.2. Nonlinear Interactions in Intermediate Turbulence 

We have deduced the form of the spectrum of intermediate MHD turbulence from 
scaling arguments based on 3-wave couplings. Moreover, these interactions are weak in 
the sense that A^a ^ 1- This might suggest that 3-wave interactions dominate those of 
higher order and that a rigorous derivation of the steady-state cascade might result from 
truncation at this order. Unfortunately, this is incorrect; interactions of all orders have 
similar strengths. 
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Consider the distortion suffered during a single collision between oppositely directly 
wave packets of similar strength having perpendicular and parallel dimensions A ^ £ and 
i. We assume that Ivx <^ Af^. Contributions from interactions involving n- waves may be 
written as 



Clearly 3-wave interactions dominate those of higher order for individual collisions. 

Next we consider the cumulative distortion due to n-wave interactions as a wave packet 
travels a distance £ <^ z ^ L*. We can picture the distortion as arising from the shearing 
of the packet as it follows the differential wandering of neighboring field lines.0 The net 
displacement of individual field lines over distance z defines a vector field whose shear 
tensor transforms the packet's shape.0 For £ <^ z <^ this transformation is close to the 
identity, so it may be expanded in a Taylor series. The expansion parameter is {z/L^Y^"^, 
the dimensionless measure of fractional spreading over distance z of a bundle of field lines 
whose cross sectional radius at 2; = is A. Terms of order n — 2 in f a correspond to n-wave 
interactions.0 These terms have the form 

S'^Vx ( (iVxY'\ Z\i^-^)/^ ^^^{n-2)/2 ^^^^ 



A notable feature of equation (|T7|) is that contributions from higher order (n > 4) 
interactions carry extra factors of {z/lY^"^. By extra factors we mean those beyond the 
single factor {z/tj^^"^ expected to arise from the addition of a sequence of independent 
interactions between pairs of wave packets. These factors reveal an interdependence among 
collisions associated with the nonlocalized field line displacements of the wave packets.^ 

Intermediate turbulence is nonperturbative because distortions of all orders become 
large as z ^ L^. This is as expected because the cascade time across scale A is tA ~ L^/va- 



uniform displacement of a wave packet does not contribute to the energy cascade. 

■^This transformation is subject to the constraints of fluid incompressibihty and magnetic flux freezing. 

^Since wave packets foflow field lines only to lowest order, nonkincmatic terms appear at orders n > 4. 
Those of n = 4 are given in §3.4. 

^These extra factors of (z/^)^/^ are absent in the corresponding formula for wave packets which are 
localized in displacement. 
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It implies that n-wave interactions of all orders n > 3 make comparable contributions to 
the energy cascade. 



A fluid element whose Lagrangian coordinate is a has Eulerian location x, at time t, 
given by 



where ^(a, t) is the displacement field. Velocity and magnetic field perturbations at the 
Eulerian location are defined in terms of the displacement vector by 



SG employed a formulation of the MHD action due to Newcomb (1962), and developed a 
Lagrangian perturbation theory for weak MHD turbulence. The Lagrangian is a functional 
of the strain tensor field, in other words, the gradient of the displacement vector field. 
Expansion of the Lagragian density in powers of the strain tensor yields terms of second 
order, n = 2, and then fourth and higher orders, n > A. The absence of third order terms 
signifies that wave packets follow field lines to lowest nonlinear order. The absence of 
Lagrangian perturbations based on 3-wave interactions implies that Eulerian perturbations 
due to 3-wave interactions are purely kinematic. Kinematic contributions to Eulerian 
perturbations are also present at each higher order, n > 4. But these are augmented by 
dynamic perturbations arising from order n > 4 terms in the expansion of the Lagrangian. 

As before, we consider small amplitude, (ivx/Xvi) -C 1 wave packets localized in v and 
b but not in ^. Roughly speaking, convergence of the perturbation expansion requires the 
components of the strain tensor to be smaller than unity.Q An ensemble of independent 
wave packets generates an energy spectrum which is fiat for kzi ^ 1. Since u = VA\kz\, 
the power spectrum of the displacement vector field varies as for kzi <^ 1. The same 
behavior characterizes the power spectra of some of the components of the strain tensor. It 



^From this point on it is best to proceed in Fourier space, since the breakdown of perturbation theory 
is closely related to the behavior of the energy spectrum (i.e. velocity, or magnetic field power spectra) at 
small kz- 



2.2.3. Lagrangian Perturbation Theory 



X = a + ^(a, t). 



(18) 





(19) 
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implies that these components diverge. The divergence is the mathematical expression of 
the spreading of field line bundles described in §2.2. In this light, the failure of perturbation 
theory is seen to be generic, and not just a consequence of an unfortunate choice of 
perturbation variable. 

Next we investigate the effect of cutting off the energy spectrum below k^L ~ 1, where 
L £. The most strongly divergent components of C,ij have power spectra given by 

Multiplying equation ( ]20|) by k'j_ and integrating from kz ~ L^^ to kz ~ we obtain the 
average value of due to power in this wave number interval, 

E7F'^A-/;;*.fc(k)p^(g'f. (21) 

Thus 

EF~^, (22) 

where is defined in equation (0). Once again we see the crucial role played by L*; the 
perturbation expansion converges if the energy spectrum is cutoff below kzL^, ~ 1. 

The absence of third order terms in the expansion of the Lagrangian signifies the 
absence of resonant 3-wave interactions in perturbative MHD turbulence. Weak MHD 
turbulence based on 4-wave interactions is discussed in 5 3.2 



3. ON WEAK AND INTERMEDIATE TURBULENCE 

3.1. Types of Turbulence 

We have discussed, at some length, the intermediate cascade. It is time to state in a 
precise manner the properties that characterize the three kinds of turbulence. Let us begin 
with a definition. Weak turbulence is characterized by the following properties:^ 



^For a standard discussion of these points see the introduction to "Kolmogorov Spectra of Turbulence P 
by Zakharov, L'vov, and Falkovich, 1992. 
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I. Nonlinear interactions among an ensemble of waves which are weak in the sense that the 
fractional change in wave amplitude during each wave period is small. 

II. The existence of a convergent perturbation expansion for the nonlinear interactions. 
Typically the small parameter is the fractional change in wave amplitude during a wave 
period. 

When these two conditions are satisfied, a formal theory of resonant wave interactions may 
be derived. For the theory to be nonempty, either the 3-wave or 4-wave resonance relations 
must possess nontrivial solutions. Power spectra of cascades arise as stationary solutions of 
the kinetic equation describing modal energy transfer. 

Specifically, for MHD turbulence, we find 

1. The turbulent cascade based on 4-wave interactions, derived in SG, is the unique weak 
cascade that satisfies both I and II above. Moreover, this cascade is realizable; it could 
in principle be set up experimentally. While 3-wave interactions do not vanish, they are 
nonresonant and do not transfer energy among different waves. 

3. The critically balanced cascade described in GS is an example of strong turbulence. 
It violates both conditions I and II. The interaction time is of order the wave period. 
Interactions of all orders have comparable strengths; there is no valid perturbation 
expansion. 

3. There is a third type of MHD turbulence that satisfies I, but not II. It is an example of 
what we have called Intermediate Turbulence. This turbulence exhibits weak interactions, 
but strains in the fluid are so strong that perturbation theory diverges. The 3-wave 
interactions are included, but not dominant; it turns out that interactions of all orders 
contribute equally weakly. The inertial range spectrum of intermediate MHD turbulence is 
given, for the first time, in equations (110) of this paper. 



3.2. Weak Alfvenic Turbulence 

Having discussed intermediate turbulence in some detail, we provide a brief outline 
of the theory of weak MHD turbulence that SG constructed using resonant 4-wave 
interactions. By studying the resonant terms of the fourth order Lagrangian, SG derived 
a formal kinetic equation for the evolution of energies (more precisely, "wave action") in 
different modes, and proved that a cascade of energy, indeed, emerges as a stationary 
solution. The elementary interactions involve scattering of two waves which respect the 
following conservation laws: 



ki + k2 = ka + k4 , cui + co'2 = c^s + 004 



(23) 



Using ujk = VA\kz\, and the z component of the equation involving the k's, SG prove that 
kiz = ^32 > , and k2z = k^z < . Of course, the symmetry of equations ( p3D allows us 
to flip the signs of all the kz, or permute indices 3 and 4; the important point is that the 
scattering process described by equations (pSf) leaves the kz components unaltered. This 



implies that waves with values of kz not present in the external stirring, cannot be created 
by resonant 4-wave interactions. 

In addition to developing a formal theory, SG also provided a heuristic derivation of 
the weak 4-wave cascade for shear Alfven waves. Here we note the main properties of the 
weak cascade of shear Alfven waves: 

(i) As discussed above, there is no transfer of energy to small spatial scales in the z direction; 
the energy cascade in k-space occurs only along k^ (i.e in directions perpendicular to the 
mean magnetic field). 

(ii) The three dimensional energy spectrum, E, is defined by 



E^I = /m,fcx)|^, (24) 

where X) is a sum over wave packets of various scales. Weak turbulence relies on a 
convergent perturbation theory. As discussed in § 2.2.3 this requires that the spectrum, E, 
be cut off for j/Cz-^*! < 1 • Moreover, weak 4-wave interactions do not change kz, which 
implies that E may have a quite arbitrary dependence on \kz\. This simply depends on the 
nature of the excitation, and is not of much interest here. If A ~ /cj^ is a perpendicular 
length scale belonging to the inertial-range, the scalings derived by SG for the weak, 4-wave 
cascade are 



(iii) The number of collisions needed for the packet to lose memory of its initial state. 
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Note that, in common with the spectrum of intermediate turbulence (given in equation 10), 
A^A decreases as the cascade proceeds to smaller A. 

(iv) In their treatment of weak turbulence, SG unwittingly made the assumption that E 
was cut off at small \kz\. While the weak turbulence of SG is realizable (see § 4 later in this 
paper), this feature makes it less applicable than intermediate turbulence. 



3.3. A Controversy, and its Resolution 

NB have proved, analytically, that 3-wave interactions between small amplitude wave 
packets are non zero if the kz = Fourier components are non zero. On the other hand, 
SG have argued, using a Lagrangian perturbation theory, that 3-wave interactions are 
absent in weak MHD turbulence. In § 2 we claim that the interactions found by NB lead 
to intermediate, rather than weak turbulence. Here we bolster that claim by an explicit 
evaluation of 3-wave and 4-wave interactions in Lagrangian coordinates. 

The simplest derivation employs the Lagrangian displacement vector field as the basic 
variable (cf. eqns. and |]19|)- Incompressibility implies that the transformation between 



Lagrangian and Eulerian coordinates have unit Jacobian. Thus 



J = 1 + V ■ ^ - i : + i ( V ■ I) V ^ce'e' - ^ (V ■ I) ce' + ^ (v ■ , m 

where V refers to derivatives with respect to Lagrangian coordinates, and 



: = e^e'- (28) 

The displacement vector ^ is split into transverse and longitudinal components, t] and 
C respectively, such that 



^ = r7 + C, (29) 

with 

V ■ T7 = . (30) 

The components of r] are the two independent variables; C is obtained from equation (p7|). 
Thus 
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(31) 



and 



(32) 



The first term on the right hand side of equation (^) is included for completeness since 
r}2 = 0.0 

Following the development in § 3 of SG, we write the Lagrangian as 



dt 



daw 



(33) 



However, we present our calculations in real space, rather than in Fourier space; the results 
are identical, but the real space version turns out to be useful later. Solving equation (|3T|) 
yields 



C2 = -v 



Stt la — a'l 



(34) 



We now write 



C — C2 + jCi, 



(35) 



(the third order terms vanish) with 



£2 = ? / d^a 



dt 



drji 

dan 



(36) 



and 



£4 = L / d-'a 



dC2 
dt 



- V, 



dC, 

daw 



(37) 



^See equation (^) and the footnote which follows it. 
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Variation of the action 

S = j dt\c + j (faP{y-r))^ (38) 
with respect to 77 leads to the (Euler-Lagrange) equation of motion.^ 

3.3.1. 3- Wave Interactions In Lagrangian Coordinates 

In this subsection we recover the results of the 3-wave interactions calculated by NB. 
Moreover, we demonstrate that they are purely kinematic in Lagrangian coordinates. 



To lowest order, the contribution of £4 may be ignored. Using equation ( p6D for C2 in 
the variation of 5* in equation results in the following simple, linear equation: Q 

whose general solution is a superposition of wavepackets traveling in the positive, and 
negative z directions: 



= 77+(a_L, ay - VAt) + rj^ {a±, ay + VAt) . (40) 

NB calculated, in Eulerian coordinates, the lowest order perturbation due to a collision 
between oppositely directed wave packets. It is a trivial matter to obtain this quantity 
by using Lagrangian perturbation theory. To do so, we transform the right side of the 
expression for v, given in equation (p!9|), into Eulerian coordinates. To first order, 

vi(x,t) = ^(x,t), (41) 

is the velocity field of unperturbed wave packets, and 77^ is the corresponding displacement 
field. To second order, 



is a Lagrange multiplier needed to insure that V • T7 = 0. 
^°We assign first order to rj of unperturbed wave trains. 
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V2 = -(r7i-V.)vi + ^, (42) 



where the subscript 'x' refers to Eulerian coordinates.|^ When we substitute equation (^ 
on the right side of equation (02|), we obtain three different types of terms; those that 



contain two powers of rjf or rji have nothing to do with perturbations induced by colhsions. 
Only the mixed terms describe distortions suffered by a wave packet during a colhsion with 
an oppositely directed wave packet. If 



Arif = T[vt{^±,+oo)-rif{x^,-oo)], (43) 

is the net displacement of field lines due to the +/— wave packets, then the asymptotic 
distortion due to a collision is 

Av2± = - Ar/f • V, v± - V,. / II! - I 44 

J An |x — x'l 

The above expression is equivalent to equations (15) and (16) of NB. The wave packet 
distortions expressed by equation (|4^) are kinematic, as described in § 2.2.3. Displacements 
resulting from any sequence of collisions are obtained by summing individual values.0 Note 
that both terms on the right side of equation (S) depend on Arj^, the net displacement of 
field lines, which is related to the Fourier amplitude of rjf with = 0. 



3.3.2. ^-VFaue Interactions in Lagrangian Coordinates 

Our principal aim in this subsection is to demonstrate that 4-wave interactions obey 
the distance scaling proposed in equation (^). Our secondary goals are to obtain the 
spectrum of weak Alfvenic turbulence from a configuration space calculation, and to recover 
the frequency changing terms discovered by NB in full MHD.|^ 



^^The absence of dri2/dt is due to the vanishing of £3. 

^^SG missed these distortions because they studied wave packet interactions in which the wave packets 
were locahzed in rj; for these Arjf = 0, and third order couphngs vanish even in Eulerian coordinates. This 
is consistent with the point made in footnote 5 of SG, that resonant couphng coefficients are independent 
of the variables used. Of course, for this to be true, one must be in a regime where perturbation theory is 
valid, which obtains only for wave packets localized in rj. 

^•^NB's discovery was made using reduced MHD. 
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We have shown that the 3-wave interactions of NB arise from kinematic perturbations 
in Eulerian coordinates. Dynamic perturbations require the interaction of at least 4 waves 
and are associated with perturbations in Lagrangian coordinates. We already possess the 
machinery necessary to derive these. When £4 is included in the action of equation (P^D, 
the variation with respect to r/i leads to Euler-Lagrange equations, which may be thought 
of as adding third order terms to the right side of equation (|39|) : 

where P3 is determined by requiring that V ■ 773 = . 

For simplicity, we evaluate the deformation suffered by a wave packet traveling in the 
positive an direction. It proves convenient to transform to coordinates 



a = a\\ — VAt, T = t. (46) 

As given in equation (|40|), the unperturbed wave packets in these coordinates have the 
forms 

r/l = 77|(a_L,a) r]^ = rj^ {a.±, a + 2vat). (47) 
The equation of motion for rj^ in the new variables, r and A = (a_L, a), reads 



where the primes indicate dummy integration variables, and the choice of superscripts 
applied to r]i on the right side of the equation is dictated by the requirement that the 
differential operators acting on the integral don't kill it; one r]f and one 773 must appear in 
the integral over d^A' because the derivatives d/dr and {d/dr — 2vAd/da) annihilate rjf 
and 77^, respectively.!^ 



^''Proof of the second of these relations requhes an integration by parts to transfer d/da acting on 
/|A — A'l to d/da' acting on 'Vrj^ : ^rji. 
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Further expansion of the right side of equation (^Hf), obtained by writing "H^ = rji + ri{ , 
results in two terms. Rather than carry the cumbersome pressure term through the 
remainder of our calculation, we discard it and replace 773 on the left side of the equation of 
motion by T/3 to remind us that its longitudinal part must be subtracted off at a later stage: 



d f d d \ ^ ( d d \ , _ d „ f d^A' Vr/^ : Vrjt 

7^ - 2^A7r- % 









(-- 


2VA 




[dr 


da J 










■V) 


(- 

[dr 



drXdr da I {dr da}'" ' dr J Atx |A - A'| 




An \A - A' 



(49) 



Although we have used coordinates, (a,r), moving with the positive wave packet, we 
emphasize that the equation of motion (eqn. |^^) remains symmetric in f\t^ and f\^: i.e., 
the action of the first term on a +(— ) wave packet is identical to that of the second term 
on a — (+) wave packet. To determine the perturbations suffered by a + wave packet, we 
set 7)3 = 7)3 on the left side. 

If the first term were the only perturbing interaction, we could peel off the operator 
{d/dr — Iv^dlda) from both sides. Integration over time would yield 

A*,i ^ - .r(,r ■ V)V / ^ ^'"l^^fi-^"' , (50) 

which describes the 4- wave interactions (cf. eqn. ||23|) that form the basis of the weak 
turbulence theory of SG. 

The third order Eulerian velocity perturbation consists of both kinematic and dynamic 
terms. An explicit expression follows from equation ([T9|): 



1 Qt} UTl uC 

V3(x, t) = -ri,r,, : V,V,^ + (771 ■ V,.) (ti, ■ V,) ^ - (ti, ■ V,) 

The final term is the sole dynamic entry. Each of the other terms is constructed from rj^ 
without use of the equation of motion. In particular, ^3 is obtained from equation ( ^2]) 



Let us estimate the fractional distortion suffered by a positive wave packet after 
traveling a distance z ^ i. Among the plethora of kinematic terms, consider only 
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Av3+ = (r;rr/r) : V.V.^ . (52) 



2 V-'^ -'W • • ^ • ^ dt 

To order of magnitude both this term and the dynamic term contribute 



(53) 



the same order of magnitude as the fractional Eulerian distortion given by equation (|TTp.Q 

To make contact with the weak 4-wave cascade, we suppose that the velocity spectrum 
of the negatively directed waves is cut off at ~ L-^ < £-^[3 Then the distort ions build 
up coherently over distance L. Over longer distances, z ^ L, they add with random phases 
implying 

Cascade occurs when 5^v\/vx ~ 1; thus 



-l) ■ <=^' 



But for the extra factor of (i/L), this is identical to equation (p6D, whereupon, provided 
A^A ^ 1, the spectrum of the weak 4-wave cascade, given in equation (pS]) follows. 



Next we return to investigate the second term in equation (|49|). Imagine that the first 
term were switched off. Then remove a d/dr from both sides leaving the left side in the 
form {d/dr — 2vAd/da)fj^i,. We expect the action of d/dr to be sub-dominant. Therefore, 

implying that the net change in the + wave packet is 



^^We draw particular attention to the linear dependence on z/£. 

^^For simplicity we consider the symmetric situation and drop the ± signs in the superscripts. 
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We have written equation { ^7\) in a form that makes exphcit the dependence of the 
perturbation on only the net displacement induced by the — wave packet. In other words, 
the perturbation induced in a wave packet is proportional to the amplitudes of the = 
Fourier components of oppositely directed wave packets: these terms were first discovered 
by NB in the limit of reduced-MHD. We may describe the interactions by the following 
kind of resonant 4-wave process, 



ki = k2 + kg + k4 , 

UJi = UJ2+ UJ3 + ^4:, (58) 

wherein one wave may be induced (by an oppositely directed wave) to split into two, or 
two waves could equally well be induced to combine into one. Manipulating the resonance 
relations leads to kiz, k2z, ^32 positive (say), while k^z approaches zero from below. Thus 
uji = UJ2 + allows for changes in the frequencies of wave packets, and harmonics, as well 
as sub-harmonics are generated by this process. These 4-wave interactions require the 
energy spectrum of the negatively directed packets to be flat near kz = 0; thus they do not 
arise in weak turbulence. Are they of importance to intermediate turbulence? To this end, 
let us make an order of magnitude estimate of the perturbation. From equation (|57|), the 
perturbation suffered by the + wave packet upon traveling a distance z ~ v^t is 

Svt-kl7]+7]+A7]^, (59) 
where Ar]^ = r]-{\±, z) - r]-{\±, 0). Using rjx ~ i{vx/vA), 

Ar]7^^{£zy/\ (60) 
Va 

Thus the fractional distortion of the positively directed wave packet, 

5i]x ( ivx \ ( iv^ \ ( zX^/"^ 



r\j ; — 



Vx m \>^vaJ \XvaJ \iJ 



(61) 
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which is smaller than the corresponding quantity in intermediate turbulence (set n = 4 
in eqn. |T^, and compare with the above eqn. |61|]) by a factor (z/^)^/^; thus harmonic 



generation is unimportant for intermediate turbulence.]^ 



4. EVOLUTION OF A WEAK PERTURBATION 



Now we address the issue of the physical relevance of the cascade proposed by SG for 
weak turbulence. We go on to show that weak turbulence, intermediate turbulence, and 
strong turbulence can be consecutive stages of a single turbulent cascade. 

Imagine exciting shear Alfven waves (isotropically on scale £/L <^ 1 with amplitude 
vi/va <^ 1) in a cubical box filled with an electrically conducting fluid which is threaded by 
an unperturbed magnetic field aligned parallel to the z axis. Let the box have side length 
L and assume that its walls are made of an excellent electrical conductor. Then field line 
displacements associated with turbulent motions must vanish at the walls. This boundary 
condition provides the cutoff of the energy spectrum for kzL < 1. 

Suppose that 



— «-. 

vaJ L 



(62) 



Then field lines initially separated by scale I maintain this approximate spacing. This 
ensures that there are no resonant 3-wave interactions in the upper part of the cascade and 



that resonant 4-wave interactions dominate. From equation (|2^) we have 



va 




(63) 



The weak cascade grades into the intermediate cascade when 



which occurs for 



(64) 



However, it might play an important role in strong turbulence. 
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For the intermediate cascade 

VA \vaJ \iJ \ij ' 

The intermediate cascade steepens into the strong cascade when 

Vx A 

VA r 

which takes place at 



e \vaJ \i. 

Within the inertial range of the strong cascade 



VA \Va) \d V) ' 

The complete three dimensional, inertial range spectrum is given by 



(65) 



(66) 



(67) 



(68) 



(69) 



' (A;x^)-^°/^ if (L-i <h< i-\ and i-^ < k± < 
(M)"^ ii{L-^<k,<i-\andX^^<k^<X^^); (70) 
{k±ey^/^, if (L-i <k,< e-\ and ki_ > A2 ^ ). 



There are a couple of points worth noting in connection with the above combined cascade. 
The intermediate cascade is confined to A2 ^ A ^ Ai where the ratio 

is independent oi vi/va- As 7?. — > 1 from below, the inertial range of this cascade shrinks 
to zero exposing a direct transition between the weak and strong cascades. This is the 
transition discussed in SG and GS. 
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5. DISCUSSION 

In a recent paper, Ng & Bhattacharjee (1996) claim that (i) in weak MHD turbulence, 
3-wave interactions between oppositely directed wave packets are non zero if the = 
components are non zero; (ii) 3-wave interactions dominate over 4-wave interactions. 
We hope to have persuaded the reader that (i) is correct so long as "weak" is altered to 
"intermediate" , and (ii) is true only for individual collisions between small amplitude wave 
packets. However, NB deserve full credit for demonstrating the importance of 3-wave 
interactions, as well as discovering the frequency-changing terms in 4-wave interactions. 
Both are a consequence of a non zero net displacement of field lines, due to the perturbations 
of wave packets that have localized velocity and magnetic field perturbations. 

Intermediate turbulence shares with weak turbulence the property that wave packets 
are long lived; interaction times are much longer than the wave periods. However, 
the distinguishing feature is that perturbation theory in not applicable to intermediate 
turbulence; this should be clear from our demonstration that interactions of all orders 
have the same strength. Of course, during individual collisions, 3-wave interactions 
dominate over all higher order interactions. However, as described in § 2.4, for wave packets 
localized in velocity and magnetic fields, but not in the net displacement of field lines, 
subsequent collisions are correlated; this makes all n-wave interactions contribute equally 
to intermediate turbulence. If Iroshnikov and Kraichnan had performed their heuristic 
estimates, taking account of the fact that there is no parallel cascade for long lived wave 
packets, they would have found the spectrum of the intermediate cascade. This is one 
case in which the assumption of isotropy (usually innocuous) is misleading; the anisotropic 
intermediate cascade strengthens on small spatial scales, whereas the isotropic IK cascade 
weakens. 

We have devoted attention almost exclusively to shear Alfven {sA for brevity) waves, 
ignoring the dynamics of pseudo Alfven {pA for brevity) waves.Q A generic excitation 
may be expected to put power equally into both kinds of waves. Shouldn't we then study 
the interactions of the pA waves among themselves, as well as with sA waves? Will this 
modify the cascades we derived for the sA waves? It turns out that the pA waves are 
slaved to the sA waves. The reason is as follows. Suppose that we were following the 
distortions suffered by a positively directed wave packet due to other, negatively directed 
ones. The nonlinear interactions are given, to order of magnitude, by some power of 
(v~ ■ V)v"*" ~ (k • v~)v+. Because the cascades in MHD turbulence are anisotropic, we 



^^The only exception is the proof that there are no 3-wave coupHngs for either type of Alfven wave, in 
the Lagrangian perturbation theory of SG. 
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expect that k± 3> kz in the inertial range. We note that the polarization of an sA wave 
is essentially along ixk^, and that of a pA wave is along z, and hence estimate that the 
"operator" k ■ v ~ {k±VsA) + (k^VpA) ~ (/i;±fs^)-Q Thus the pA waves are slaved to 
the sA waves; in a later paper, we will derive kinetic equations, and demonstrate that the 
spectrum of the pA waves is identical to that of the sA waves for weak, intermediate, and 
strong turbulence. 

How might an energy spectrum with a cutoff below fc^L ~ 1 arise? Two related 
possibilities come to mind. The first involves a thought experiment that could be realized 
computationally; we described this in § 4. A more natural setting might be the atmosphere 
of a massive star which possesses a strong external dipole field.Q The energy spectrum of 
waves in the stellar magnetosphere would be cutoff on scales longer than the length of the 
flux tubes that link the northern and southern magnetic hemispheres. The cut off at small 
kz necessary for weak turbulence makes it, in general, less applicable than intermediate 
turbulence. However, the limited inertial ranges of both weak and intermediate cascades 
suggests that neither is likely to find much application in nature. The critically balanced 
cascade, proposed by Goldreich & Sridhar (1995) for strong MHD turbulence, remains as 
the most likely candidate for interstellar turbulence. 
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